B® — B® Mixing Parameters from Unquenched Lattice QCD 
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We determine hadronic matrix elements relevant for the mass and width differences, AM S & AF S , 
in the B a —B a meson system using fully unquenched lattice QCD. We employ the MILC collaboration 
gauge configurations that include u, d and s sea quarks using the improved staggered quark (AsqTad) 
action and a highly improved gluon action. We implement the valence s quark also with the AsqTad 
action and use NonRelativistic QCD for the valence 6 quark. For the nonperturbative QCD input 

into the Standard Model expression for AM S we find f Bs \J Bb s = 0.281(21)GeV. Results for four- 
fermion operator matrix elements entering Standard Model formulas for AF S are also presented. 

PACS numbers: 12.38.Gc, 13.20.Fc, 13.20.Ho 
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INTRODUCTION 

Recent developments at the Tevatron Run II have dra- 
matically improved our knowledge of the mass difference 
AM S between the "heavy" and "light" mass eigenstates 
in the B® — B® system. The Spring of 2006 witnessed first 
the two-sided bound on AM S by the D0 collaboration [l| 
followed quickly by a precise measurement of this quan- 
tity by the CDF collaboration B s mixing occurs in 
the Standard Model through box diagrams with two W- 
boson exchanges. These diagrams can be reexpressed in 
terms of an effective Hamiltonian involving four-fermion 
operators. In order to compare the Tevatron measure- 
ments with Standard Model predictions, matrix elements 
of the four-fermion operators between the B® and B® 
states must be computed. Only then can one test for 
consistency between experiment and the Standard Model 
and, in the case that precise agreement fails to be real- 
ized, hope to discover hints of new physics. B° s - B° s 
mixing is a AB = 2 process and sensitive to effects of 
physics beyond the Standard Model. Hence a large effort 
is underway to nail down the Standard Model predic- 
tions as accurately as possible. In the current article we 
present a fully unquenched lattice QCD determination 
of the hadronic matrix elements of several crucial four- 
fermion operators. 



SIMULATION DETAILS 



m s being the physical strange quark mass. For the 
strange valence quark we use the improved staggered 
(AsqTad) Q quark action. The 6-quark is simulated 
using the same improved nonrelativistic (NRQCD) ac- 
tion employed in recent studies of the T system H and 
for calculations of the B and B s meson decay constants 
0, 0] and the B — > n, Iv semileptonic form factors 
As in our previous work using the MILC configurations 
we use the T 2S-1S splitting to fix the lattice spacing, 
which in the present case gives oT 1 = 1.596(30)GeV 
and a- 1 = 1.605(29)GeV Q for the m f /m s = 0.25 
and mf/m s = 0.5 ensembles respectively. The bare b 
and s quark masses have likewise been fixed already in 
previous simulations of the T Q and kaon [§] systems. 
Some theoretical issues remain having to do with the 
need to take a fourth root of the AsqTad action deter- 
minant while creating the MILC unquenched configura- 
tions. This procedure is the focus of intense scrutiny by 
the lattice community and there has been considerable 
progress in our understanding of the issues involved dur- 
ing the past year [l(| . To date no obstacles have been 
found to invalidate obtaining true QCD in the contin- 
uum limit. The MILC configurations and the light and 
heavy quark actions employed in this article have also 
been tested by comparing the results of accurate calcu- 
lations for a large range of hadronic quantities to experi- 



mental results [ill . [12] , Il3l . Il4l [15j . The outcome of these 
tests have been very encouraging. Here we apply the 
same successful lattice approach to B® — B® mixing. 



Our simulations use the MILC collaboration Nf = 2+1 
unquenched gauge configurations To date we have 
completed calculations on two of the MILC coarse en- 
sembles with the light sea quark mass mf satisfying 
rrif/m s — 0.25 and m//m s = 0.5 respectively and with 



THE FOUR-FERMION OPERATORS AND 
MATCHING 

We have studied the following four-fermion operators 
that enter into calculations of AM S and Ar s in the Stan- 
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dard Model ("i" and "j" are color indices) 



OL = [¥ s^v-aIV S>} V -A, 



which are of the form, 



OS = [b l s^s-piV s j }s-p, 
03 ee \t? s^s-pIIp s^s-p. 



(1) 
(2) 
(3) 



In continuum QCD in the MS scheme matrix elements of 
these operators are parametrized in terms of the B s me- 
son decay constant fs 3 and so-called "bag" parameters 
B(jj) at some scale 



(B 8 \OL\B s )ffi 



\fl s B B MMl a . (4) 



The factor | is inserted so that Bb, = 1 corresponds 
to the "vacuum saturation" approximation. The four- 
fcrmion operators OS and 03 have similarly each their 
own bag parameter 



(OS)^ s 



2 B S (IJ) 2 



/ni\ MS — 1 f 2 B s(^ m 2 



with 



1 



Mr 



(m b 



(5) 



(6) 



(7) 



The Standard Model expression for the mass difference 
AM S is given by 0, 



AM, 



^|y*^| 2 ^So(, t )M Ba 44, 



(8) 



where Xt = mf/M^, r/f is a perturbative QCD cor- 
rection factor, So(xt) the Inami-Lim function and Vt s 
and Vfh the appropriate Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements. The nonperturbative QCD in- 
put into this formula is the combination f B Bb s with 
Bg a the renormalization group invariant bag parameter. 
At two-loops and using nf = 5 and a^^°\fi = m b — 



4.8GeV) = 0.212 [151 one finds B B JB Bs = 1.534. 

In order to evaluate hadronic matrix elements of the 
four-fermion operators via lattice QCD methods, one 
must first relate the operators in continuum QCD to op- 
erators written in terms of lattice heavy and light quark 
fields. We carry out this matching between continuum 
QCD and the lattice theory through 0(a s ), 0(^p-) 
and O(^). Our lattice theory works with NRQCD b- 
quarks. At lowest order in 1/M the b fields in - ([3]) 
must be replaced by NRQCD heavy quark or heavy anti- 
quark fields. The tree-level relation between NRQCD 
and full QCD fields is given bv^ the Foldy-Wouthuysen- 
Tani transformation. At 0( this brings in di- 

mension seven corrections to the four-fermion operators, 



OLjl 



M-{^ 

+ [6*sV-a[V&j-7s j V-a}. ( 9 ) 



Similar 1/M corrections OSjl and 03 jl can be intro- 
duced for the four-fermion operators OS and 03. To 
the order stated above, matching between (OX) MS (X 
= L,S or 3) and matrix elements in the lattice theory is 
then given by (we suppress the fi dependence), 



2M f 



(OX) 



MS 



[l + a s - p XX ](OX) + a s ■ Px y(OY) + 
[(OXjl) - a s ({* x (OX) + a Y (OY))} . (10) 

(OX) without the superscript MS stands for the matrix 
element in the lattice theory. Even at lowest order in 
1/M there is mixing between the four-fermion operators. 
At 0(a s ) the mixing occurs between X, Y = L and S for 
(OL) and (OS) and between X, Y = 3 and L for (03). 
This mixing takes place already in continuum QCD when 
one carries out an expansion in 1 /M [13, EH ■ Due to the 
good chiral properties of AsqTad light quarks, which we 
use for the valence s quark in our simulations, no addi- 
tional operator mixing arises upon going from the contin- 
uum to the lattice theory. We have multiplied (OX) MS 
in (p~0|) by a factor of 2 £ { in order to take into ac- 
count the different normalization of states in QCD and 
the lattice theory and also to render the lattice matrix 
elements (OX) dimensionless. Details of calculations of 
the one-loop coefficients pxY and ( X Q Y will be presented 
in a separate paper. The methodology is similar to that 
of [11, [H. As in those matching calculations for heavy- 
light currents, the a s ■ Qi x and a s ■ Cio Y terms in (fT0|) 
are necessary to subtract O(-^jj) power law contributions 
from the matrix elements (OXjl). 



SIMULATION RESULTS AND ERROR 
ESTIMATES 

The hadronic matrix elements (O), 0=OX and OXjl, 
are determined by evaluating three-point correlators via 
numerical simulations, 

c^\hM) = E (0|*5.(iMi) [d](o) -t 2 )|o). 

(11) 

<f>B s is an interpolating operator for the B s meson and 
the four-fermion operator O is fixed at the origin of the 
lattice. We fit C^ 4 ^ together with the B e meson two- 
point correlator, C B (t), to the following forms 
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TABLE I: Matrix elements in the lattice theory for fixed 6 and 
s valence masses and two values of the light u, d sea quark 
mass. Errors are statistical plus fitting errors. 



TABLE II: Error budget for quantities listed in (TPo)) 





rrif /m s = 0.25 


m//m s = 0.50 


(OL) 


0.1036(83) 


0.1069(92) 


(OS) 


-0.0680(54) 


-0.0687(61) 


<03) 


0.0143(12) 


0.0142(13) 


(OLjl) 


-0.0227(18) 


-0.0229(18) 


(OSjl) 


-0.0130(10) 


-0.0139(11) 


(03jl) 


0.0021(3) 


0.0026(3) 


(OLjl) (aub) 


-0.0138(20) 


-0.0140(20) 


(OSjl) {sub) 


-0.0072(11) 


-0.0081(12) 


(OSjl) {sub) 


0.0008(3) 


0.0012(4) 



N cxp -1 
j,fc=0 



(12) 



C B {t) = X)<0|* B .(x,t)<(0)|0> 



3=0 



j-t e -E^>{t-l)_ 



(13) 



The dimensionless matrix elements entering the right- 
hand side (RHS) of (fT0|) are then given by, 



(O) 



Ago 



(14) 



Results for (O) are summarized in Table I for our two 
dynamical ensembles. The errors are combined statistical 
and fitting uncertainty errors. More details on our fits 
In Table 1 we also show results for 



21|. 



are given m 

(OXjl)( sub ), the true relativistic corrections (after power 
law subtractions) from the dimension seven operators. 
By considering (OXjl) t SU b) / (OX) , one finds the physical 
0(A QCD /M) contribution to be -13% for (OL), 11% for 
(OS) and 6-8% for (03). 

Having determined the matrix elements in the lattice 
theory we can plug the numbers into the RHS of (fTT))) . 
For this matching we use a s — a < "y , ~ 3 \2/a) = 0.32 [lit . 
We set the scale for a s to q* = 2 /a, which is close to 
q*'s evaluated for heavy-light currents using other heavy 
and light quark actions. The matching coefficients pxy 
are generally functions of the MS scale /i through the 
combination log(-^-). We present results for fi = m b . 

We evaluate the RHS of ^ for (OL) m , (OS) m and 
(03) MS and combine with the definitions in (0| - © to 
obtain the main results of this article, namely, 
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Statistical + Fitting 


9 


% 


Higher Order Matching 


9 


% 


Discretization 


4 


% 


Relativistic 


3 


% 


Scale (a~ 3 ) 


5 


% 


Total 


15 


% 



The main errors in these quantities are listed in Table 
II. One sees that the two dominant errors are due to 
statistics + fitting and higher order matching uncertain- 
ties. We have also included a nonnegligible error coming 
from the uncertainty in the scale (lattice spacing) for the 
MILC ensembles used. At the final stage of extracting 
results for (TPo]) . one has to convert a 3 /^ s Mb s into physi- 
cal units. An uncertainty of — 1.8% in the lattice spacing 
turns at this point into a ~ 5% uncertainty for a -3 . The 
leading discretization error in the actions employed here 
comes in at 0(a 2 a s ) ~ 2% and is believed to be domi- 
nated by taste-changing effects in the Asqtad action, an 
assumption that has been checked recently by comparing 
Asqtad valence quarks with more highly improved stag- 
gered quarks from the HISQ action [22|. We multiply the 
2% by a factor of 2 to come up with a total discretization 
uncertainty of 4%, which should cover taste-changing ef- 
fects in the sea as well, including the fourth root. This 
total error is consistent with scaling tests carried out via 
explicit simulations at two lattice spacings of other B 
physics quantities such as decay constants and semilep- 
tonic form factors employing the same actions as in the 
present article 6, 8]. In the future we plan to repeat the 
current calculations at finer lattice spacings in order to 
reduce discretization uncertainties and also to carry out 
tests with HISQ instead of AsqTad light valence quarks. 

In Table II we take the operator matching error to be 
1 x a 2 since matching is done directly for the combination 
/|> s Bb b (and for the other quantities in (fl5|) ). A naive at- 
tempt to deal separately with f B and Bb, in the formula 
for AM S could increase the error estimate since the per- 
turbative error for just fs s alone (unsquared) is usually 
also taken as 1 x a 2 coming from higher order matching 
of the heavy- light current. We avoid unnecessarily sepa- 
rating out the bag parameters and possibly introducing 
ambiguities in error estimates by always working with 
the relevant combination f B Bb s ■ Several years ago ref- 
erence t 23] also emphasized the virtues of working with 
physical combinations and never splitting off the bag pa- 
rameters. The possibility of reducing errors by focusing 
on the combined f B Bs q is also mentioned in [24j |. 

Table III gives our final values for the square root of the 
quantities listed in (|15p together with the scale invariant 



(15) combination fs B y Bb s ■ One sees that the light sea quark 



4 



TABLE III: Results for the square root of quantities listed in 
(|15|) . The unhatted bag parameters are given at scale n = m&. 
Errors quoted are combined statistical and systematic errors. 
Note that percentage errors here are smaller than those given 
in Table II by a factor of two due to the square root. 



/s a v Bb s [GeV] 
fB.y/B B ,(m b ) [GeV] 



h 



•J B s (m b ) 



[GeV] 



fB y Bs ( mb) [GeV] 



mf/m s — 0.25 m//m s = 0.50 



0.281(21) 
0.227(17) 



0.289(22) 
0.233(17) 



0.295(22) 



0.305(23) 



0.301(23) 



0.310(23) 



mass dependence is small and not statistically significant. 
The largest difference between the central values of the 
mf/m s = 0.25 and mf/m s = 0.5 results is less than 3%, 
smaller than any of the other errors, and in particular sig- 
nificantly smaller than our current statistical errors. Any 
reasonable estimate of chiral extrapolation (in m sea ) un- 
certainties will not affect the total error in Table II. In 
the future we plan to use Staggered Chiral perturbation 
theory (SChPT) [25|, [26| to extrapolate to the physical 
chiral limit. SChPT formulas for B q mixing with Asqtad 
light quarks are being worked out by Laiho and Van de 
Water (2t| and will be important in Bd mixing studies 
where one needs to extrapolate in both the valence and 
sea light quark masses. In the present case of B s mixing 
and until our statistical errors have been further reduced 
and more data points are available, we do not believe the 
whole machinery of SChPT is crucial. We note that in 
our studies of the B s meson decay constant, where we 
have data for fg s at four different light sea quark masses 
on coarse MILC ensembles and for two sea quark masses 
on the fine MILC lattices, no sea quark mass dependence 
was observed @, [28[ . We do not attempt a chiral extrap- 
olation in the light sea quark mass with the current data 
and take the mf/m s — 0.25 numbers as our best deter- 
minations of the hadronic matrix elements. In particular, 

this gives for the combination fg a \J Bb s , the crucial non- 
perturbative ingredient for AA/ S , the value quoted in the 
abstract: 



1b. V B Bs = 0.281(21) GeV. 



(16) 



RESULTS FOR THE MASS DIFFERENCE AM S 

AND FOR \VZV a \ 

Using (JTHJ) one can now attempt a theory prediction 
for AM S based on the Standard Model. We plug in 



standard values for the other ingredients in ([8]) taken 
from recent reviews. We use — 0.551(7), m t (m t ) — 
162.3(2. 2)GeV (which leads to S (x t ) = 2.29(5)) and 
\V t *V tb \ = 4.1(1) x 10" 2 [H [U together with fjjfy to 
obtain 



AM s (theory) = 20.3(3.0)(0.8)ps~ 



(17) 



The first error is the 15% error from 1 B Bb s and the 
second an estimate of the error from |V t *Vtb| and m t . 
AM s (theory) is consistent with the CDF measurement 
of 17.77 ±0.10 ± 0.07PS" 1 0. At the moment the- 
ory errors, in other words the lattice errors, dominate. 
However, one is already in a position to place nontrivial 
constraints on beyond the Standard Model effects. An al- 
ternative way to test consistency of the Standard Model 
is to use the CDF measurement of AM S together with 
(fl"6j) to determine \V t * s V tb \. One finds 



\V*V tb \ = 3.9(3) x 10- 



(18) 



consistent with the value used above which is based on 
the measured value for \V c b\ phis unitarity [30(. The error 

in (|18p comes entirely from the uncertainty in fg 3 



RESULTS FOR THE WIDTH DIFFERENCE Ar s 



We have emphasized our result ([Po|) since this enters 
into AM S for which a precision experimental measure- 
ment exists. The other entries in Table III for the four- 
fcrmion operators "OS" and "03" are relevant for the 
width difference Ar s 31, for which experimental 
errors are currently still greater than 50%. A recent 
measurement by the D0 Collaboration [33[ gives, for in- 
stance, Ar s = JQ.13 ± O.O^ps- 1 . On the theory side, 
the authors of [32J have recently shown that by going 
to a new basis employing operators "OL" & "03" , as 
opposed to the old basis of "OL" & "OS", theoretical 
uncertainties from 1 /nib and a s corrections can be signif- 
icantly reduced. We insert our results from Table III for 

Jb s \J Bb s and 1b s ^-j^, takin g ag ain the mf/m s = 0.25 
data, into eq.(51) of reference [321 ] and obtain, 



Ar s = 0.10(3)ps" 



(19) 



which is consistent with the D0 measurement. About 
half the error in (|19|) comes from lattice errors in f B Bb, 

and f B ■§§■ and the other half from remaining theoretical 



uncertainties in the formula of [32[ . Both types of errors 
can be reduced through further work on the lattice and 
we plan to focus on obtaining a more accurate Standard 
Model theory prediction for Ar s while waiting for the 
experimental measurements to improve as well. 

We argued above that there is no need to separate 
out the bag parameters when calculating AM S or Ar s . 
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One might nevertheless be interested in doing so to 
judge how close we are to the "vacuum saturation" ap- 
proximation. To convert to bag parameters we use the 
central value of the B s meson decay constant deter- 
mined in 0], f Bs = 0.260(29)GeV, together with (for 
l/R 2 ) H m b = 4.25GeV and m s = 85MeV. One finds 
B B3 (m b ) = 0.76(11), B s (m b ) = 0.84(13) and B s im h 
0.90(14), values consistent with earlier quenched 1 
and Nf = 2 j35| lattice determinations. 



SUMMARY 

This article presents full QCD results for hadronic 
matrix elements of four-fermion operators relevant for 
B® — B® mixing. We give unquenched results includ- 
ing the effect of 2 + 1 flavors of sea quarks and calcu- 
late both leading and next-to-leading matrix elements 
in a nonrelativistic expansion of four-fermion operators 
for the first time. Using our nonperturbative QCD re- 
sults one finds agreement between the recent Tevatron 
measurements of the mass and width differences AM S & 
Ar s and Standard Model predictions. Our dominant er- 
rors for f B Bb 3 come from statistics + fitting and from 
higher order operator matching uncertainties. Work is 
underway aimed at reducing these errors. This will al- 
low even tighter constraints on any beyond the Standard 
Model effects entering B® — B® mixing phenomena. We 
have also initiated studies of B^ — B% mixing. Our goal 
there is to obtain precision results for the important ratio 
Ib Bb s I fB d ^B d , which will provide further consistency 
checks on the Standard Model and give us a handle on the 
CKM matrix element |Vtd|. Some errors listed in Table 
II will cancel almost completely in this ratio, such as the 
a -3 and the higher order matching uncertainties. Oth- 
ers, such as the statistical and fitting error, will cancel 
partially. 
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